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We calculate suppression of inter- and intralayer superconducting currents due to equilibrium phase 
fluctuations and find that, in contrast to a recent prediction, the effect of thermal fluctuations cannot 
account for linear temperature dependence of the superfluid density in high-T c superconductors at 
low temperatures. Quantum fluctuations are found to dominate over thermal fluctuations at low 
temperatures due to hardening of their spectrum caused by the Josephson plasma resonance. Near 
T c sizeable thermal fluctuations are found to suppress the critical current in the stack direction 
stronger, than in the direction along the layers. Fluctuations of quasiparticle branch imbalance 
make the spectral density of voltage fluctuations at small frequencies non zero, in contrast to what 
may be expected from a naive interpretation of Nyquist formula. 
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One of important problems in layered high-T c cuprate 
superconductors that is still under discussion is the origin 
of the observed linear low-temperature dependence of the 
superfluid density (phase stiffness). The latter is directly 
related to magnetic penetration depths, Ay and Aj_, for 
a magnetic field screened by currents flowing in direc- 
tions parallel and perpendicular to the superconducting 
planes, respectively Q]. This dependence is usually at- 
tributed to contribution of quasiparticles near the nodes 
of the d-wave gap. According to the alternative expla- 
nation suggested in refs. |^],[| the linear decrease in the 
temperature dependence of 1/A 2 is induced entirely by 
classical thermal phase fluctuations. Recently the role 
of fluctuations was reconsidered [Q for a d-wave super- 
conductor by means of a microscopic approach within a 
functional integral framework, and quantum phase fluc- 
tuations were found to lead to a sizeable renormalization 
of the superfluid density, the effect of thermal fluctua- 
tions being small for T <T C . 

Thermal and quantum fluctuations considered in 
refs. (|-[§ are equilibrium fluctuations, therefore, the 
problem can be solved by means of a phenomenologi- 
cal approach based on the fluctuation-dissipation theo- 
rem. Such an approach is more simple and physically 
transparent, being in the same time more general be- 
cause it is not restricted to a specific model of high-T c 
superconductivity, fn the present work we reconsider the 
role of equilibrium phase fluctuations in layered super- 
conductors applying the fluctuation-dissipation theorem 
to equations of the linear response, and come to con- 
clusions different from those of refs. Pfl- Our calcula- 
tions demonstrate that, in agreement with the results of 
rcf. at low temperatures quantum fluctuations dom- 
inate over thermal fluctuations due to Coulomb effects. 
This happens because non-uniform phase fluctuations in- 
duce a fluctuating electric field resulting in a finite energy 
of such fluctuations which is reflected in a finite value of 



the Josephson plasma frequency. So unlike refs. ^,|) we 
assert that a contribution of thermal fluctuations cannot 
account for temperature dependence of the penetration 
depth. However, in contrast to the results of ref. we 
find that renormalization of the superfluid density due to 
fluctuations at low temperatures is not large. Near T c , 
when quasiparticle density dominates over the superfluid 
density, we find a different picture. In this case ther- 
mally induced fluctuations of the quasiparticle branch 
imbalance are found to be important, and a reduction 
of the superfluid density may become sizeable. We find 
that a suppression of the superconducting critical cur- 
rent due to thermal fluctuations is larger in the stack 
direction, than in the direction parallel to the conduct- 
ing layers. This points out a possibility of a destruction 
of the superconductivity in the perpendicular direction 
at temperatures lower, than the critical temperature for 
the in-plane direction. 

The fluctuation-dissipation theorem relates correlation 
functions of fluctuations to dissipative parts of kinetic co- 
efficients describing the linear response of a system. So 
we start with the linear response of a layered supercon- 
ductor. We consider expressions for current and charge 
density written in a form of functions of gauge invariant 
vector and scalar potentials. These potentials are 

P s = (l/2)V X -(l/c)A, M = (l/2)9 tX + $, (1) 

where A is the vector potential, \ is the phase of the 
order parameter, and $ is the electric potential (we use 
units with e = 1, % = 1, and fc^ = f). The gauge invari- 
ant vector potential plays a role of the superconducting 
momentum, while the gauge invariant scalar potential is 
responsible for branch imbalance || and charging effects. 
The electric field is expressed in terms of the potentials 
as 

E = d t P s - V/i. (2) 
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In layered superconductors the variables above are re- 
lated to specific layers and must be supplied by layer 
numbers, e. g., the parallel component of P s will be de- 
noted as P n . Furthermore, the role of the perpendicular 
component of P s plays the gauge invariant phase differ- 
ence between the layers, ip n — Xn+i — Xn — 2sA±/c, where 
s is the lattice period in the stack direction. So, we may 
denote P± n = ip n /2s. 

Equations determining current and charge densities 
read 

c 2 

j||n = 4^2 P " + a 0\\dtPn - eri||V||/i„, (3) 

j±n = jc sin if n + -^-d t fn - CTl±V n p n , (4) 
2s 

K 2 

dtPn = (ld t + v h )—p n - (a 2 \\ Vjj + (J2±y 2 n )Hn + 
d t (a n V||P„ + cri_LV n (^„/2s), (5) 

where V n /i„ = (fJ> n +~L ~ M«)/ s i s the discrete version of 
the spatial derivative in the transverse direction, kT 1 is 
the Thomas-Fermi screening radius, and Vf, is the branch 
imbalance relaxation rate. 

The first terms in equations (||-^|) describe the super- 
conducting current for relative directions, and the rest 
terms are related for quasiparticle contributions. Equa- 
tion (^J) for charge density can be interpreted as a con- 
tinuity equation for quasiparticles. It was derived for 
conventional superconductors in ref. |6| and generalized 
later for layered superconductors with s-wave and d- 
wave H pairing. 

Note that the quasiparticle current does not have the 
form j = aE valid for a normal state. Instead, the re- 
sponse of a superconductor to the components of the field 
E expressed in terms of the temporal derivative of P s and 
of the spatial derivative of p (c/. (||)) are described by 
different generalized conductivities cr^, i = 0, 1, 2. The re- 
spective parts of the electric field E are induces by time 
dependent perturbations of current density, and by per- 
turbations of charge density. 

Explicit expressions for the coefficients in equations 
(HQ) depend on temperature, on scattering, and on a 
mechanism of the superconductivity (see , note that 
the explicit expressions for coefficients in ref. |^| were 
presented in the dynamic limit u> > Dq 2 where D is a 
diffusion coefficient in quasiclassic approximation). At 
temperatures T — > all conductivities are small (e. g. in 
a s-wave superconductor they are exponentially small), 
and 7 — ► 1. At T — > T c the conductivities <7j approach 
the normal state conductivities. Furthermore, the differ- 
ences between Oi for a given direction vanishes as A/T, 
7 — > A/T, h>b —>■ A/T, so that equation (||) approaches 
the continuity equation as A — > 0. Strictly speaking 
products in (Q]^) assume the convolution with respect to 
time and coordinates, i. e. in the Fourier transformed 
form the coefficients in the equations are frequency and 



wave vector dependent. 

In order to derive equations relating fluctuations Sp, 
Sip and <5P with fluctuations of charge and current den- 
sities, dp and Sj, we insert expression (Q|]) into Maxwell 
equation 

4?r 1 
VxH = -j + -<9 t D 

c c 

and into the Poisson equation. Then making Fourier 
transformation with respect to time, to in-layer coordi- 
nates r II , and to layer numbers n, we find 



(6) 
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Here q± = (2/s) sin (q±s/2), where \qj_\ < tt/s is the 
wave number obtained from the discrete Fourier transfor- 
mation with respect to layer numbers, uji — uj+iuji r , Qi — 
LU + iCli ri LUi r = AitUi^e and 0,> = 47r<7j|| are dielectric re- 
laxation frequencies. Furthermore, Co 2 — uj 2 (1 + Aj_g|), 

SI 2 = f2p(l + Ap^), where Q, p = c/X and uj p = c/Xx^fe 
are the plasma frequencies for directions parallel and per- 
pendicular to the layers, e is a dielectric constant in trans- 
verse direction, Aj_ = c/y/8wsj c . The in- layer plasma 
frequency £l p , is much larger than typical frequencies of 
the problem which are of order of the Josephson plasma 
frequency, uj p . 

Note that matrix A satisfy the Onsager symmetry rela- 
tions. Now we apply the fluctuation-dissipation theorem 
(c/. ||) to equations (||). 

From the expression for energy dissipation density Q 
which in superconductors is given by J1Q| 



Q = pd t p+jd t P s 



(7) 



one can identify potentials p and P s with generalized 
forces related to variables p and j, respectively. So ac- 
cording to the fluctuation-dissipation theorem correlation 
functions of Sjx, Sji\ and Sp are determined by imaginary 
parts of related coefficients of the matrix A in (^|). For 
example, we find 



(<y ± (q, W )<5u(q>')) = (2^) 4 (jiM(a; W)<5(q- 



(8) 



with (j 2 )^ = 2Ta' 0± /s, where T = (oj/2) coth (w/2T), 
a' = Ker. Similarly, (jxp)u = 2Tq ± a' 1± / sw, (/r% = 
2T(vi,/4ir + q^v^x + q 2(J 2\\)/ SU}2 an d so on. An alterna- 
tive way to calculate correlation functions based on the 
Langevin sources in equations of motion, and on the dis- 
sipation function of the system gives similar results. 
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Correlation functions of SP S and 6fi are related to the 
inverse matrix A~ 1 . For example, using the upper diag- 
onal component of Ar 1 we can find the spectral density 
of phase difference fluctuations 



(V)u 



^5 



j7rf(^ -^)[( 7 -|y + 
D 



(9) 



where D is the determinant of the matrix A. 

Zeros of D determine collective modes and penetra- 
tion of electric and magnetic fields into a superconductor. 
They are important in calculation of spectral density of 
voltage noise. The voltage between contacts separated 
by N layers, can be expressed according to eq.(|^) as 

N-l 

V = dt X) fn/2 — (hn ~ Mo)- Then we make Fourier 

71=0 

transformation and calculate mean square value of the 
voltage using (^J) and similar expressions for other cor- 
relation functions. We assume that the width of the 
contacts is large enough, therefore, we need functions at 
qu = 0. Finally, after some algebra we find the expression 
for the spectral density of voltage fluctuations 
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dq A 



sin 



2 qjj^ 



uj(iv-f + iv b )n 2 + eq\ [u)(lOq + lo 2 — 2u)\) — u 2 ] 
(lo 2 - louj ) (ujj+iv b )n 2 + eq\ [lo (lo 2 - lo uj 2 ) +^2^p] 
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with L = Ns. 

At low temperatures, T <C A, the relaxation frequen- 
cies are small in comparison to plasma frequencies lo p , 
and zero of the denominator gives the underdamped 
Josephson plasma mode Q|ll],[l2| at — 0. In the op- 
posite limit T> A, near T c , the plasma frequencies are 
small, flp,Ljp oc N a oc (A/T) 2 where N s is a fraction 
of the condensed electrons. On the other hand, all con- 
ductivities (Ti|| and &i± at T — » T c approach the normal 
state conductivities for parallel and perpendicular direc- 
tions, respectively. So near T c the dielectric relaxation 
frequencies are larger, than plasma frequencies for respec- 
tive directions. Furthermore, calculations for s- and d- 
wave pairing give 7 = 7rA/4T || and 7 = nA /2T^/ii7/Lu 
Q , respectively. The branch imbalance relaxation rate is 
determined by energy relaxation for s-wave pairing and 
by elastic scattering for d-wave pairing, and in both cases 
Vb oc A/T. Then zeros of the denominator in equation 
( |To| ) give the Carlson-Goldman mode |H| for direction 
perpendicular to the layers, that is underdamped in the 
case of isotropic pairing in a narrow frequency region, 
Vb,L0 2 /Lo r <C to <C (T/ A)ujp/u> r . The spectrum of the 
mode is lo 2 ps eq^j/n 2 . This spectrum in isotropic su- 
perconductors was found in clean limit in ref. Jl4| and in 
dirty limit in ref. Evolution of the spectrum from 

the Josephson plasma mode to the anisotropic Carlson- 
Goldman mode in layered superconductors was studied 
in ref. M . In the case of d-wave pairing 7 is not real, and 



the Carlson-Goldman mode is never underdamped (cf. 
§)■ 

For N ^> 1 the leading contribution to the voltage 
fluctuations is 
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This expression corresponds to the Nyquist formula and 
exhibits the Josephson plasma resonance. At zero fre- 
quency this contribution vanishes, but there is an addi- 
tional contribution to voltage fluctuations which remains 
finite at low frequencies lu <C Vb- The latter contribution 
is related to the quasiparticle branch imbalance fluctua- 
tions, and is especially pronounced near T c : 



(SV 2 ) 
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where l 2 E = Ajra^/iSbK 2 (we assumed a thick sample, 
L Ie)- Equation ( fl2"| ) contradicts to a naive interpreta- 
tion of the Nyquist theorem, according to which voltage 
noise at zero frequency is absent because the static resis- 
tivity of a superconductor is equal to zero. The non-zero 
contribution is related to a voltage drop near the super- 
conductor boundary due to penetration of the electric 
field into layered superconductor over a distance needed 
for the branch imbalance relaxation. Though (SV 2 )^ in 
( |l2| ) is independent on the distance L between the con- 
tacts, it may give a sizeable contribution which can be 
measured easily in small mesa structures. 

Using equation (|^) we can calculate mean square fluc- 
tuation of the phase difference. 



(Sep 2 ) = / (V) 



cLudq 

w 



(13) 



Since (j c sin (ip + Sip)) w j c (l — (S<p 2 }/2) sirup, the latter 
term determines renormalization of the superfluid den- 
sity. We calculate this renormalization, first, in the limit 
of low temperature. 

At T <C A the relaxation frequencies are small in com- 
parison to the related plasma frequencies lu p and fl p . We 
assume a simplifying condition ks/e> 1 which holds in 
high-T c superconductors, therefore, perturbations of Sfi 
are small and can be neglected. Then the leading terms 
in the denominator of the matrix A are 

D = U)K 2 Cll{e(ul - LU 2 )(l + Xjql) + c 2 q\ 



-4iriLj[<7 0± (l + X 2 q 2 ± ) + <7 \\e{u>l 



u; 2 + c 2 qie)/n 2 p ]}. 



A phase volume near zeros of D related to the Joseph- 
son plasma mode contributes much to the integral ( [13] ) 
resulting in a finite fluctuations at T — 0. Furthermore, 
the integral over q\\ diverges at large gii , and we cut-off it 
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at qu ~ 1/Coj where £o is the superconducting correlation 
length. In dimensional units we obtain 



To 



(14) 



3^2p s 



This result differs from that obtained in refs. For 
parameter values typical for layered cuprates we find 
To ^> T c . The magnitude of (Sip 2 ) is not large because 
the large factor A/£o in (0) is multiplied by the small 
fine structure constant. Thus the renormalization of the 
penetration depths due to phase fluctuations at low tem- 
peratures is practically temperature independent, and is 
not large. 

Calculation of (<5P||) 2 Co which determines suppression 
of the superfluid density in the in-plane direction gives 
value similar to (pr[). 

Now we calculate fluctuations at high temperatures, 
T — > T c . Since at such temperatures fluctuations are 
quasi-stationary they can be found using the standard 
approach based on functional integration of the free en- 
ergy of the system, which includes energy of the mag- 
netic field and of the superconducting current. We calcu- 
late (Sip 2 ) and (<JPi?)£o, which determine the renormal- 
ization of the superfluid stiffness in perpendicular and 
parallel directions, respectively, using again the cut-off 
at qu ~ l/£o- Then in the limit £o <C \j, which defi- 
nitely holds for layered high-T c superconductors, we find 
in dimensional units 



(Sip 2 
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Since in layered high-T c superconductors Aj^ is by few 
orders of magnitude larger, than the correlation length 
£o, the logarithm in equation ( plf ) is large as well. Then 
suppression of the critical current by fluctuations in the 
direction perpendicular to the layers is larger, than in the 
parallel direction. Using parameters typical for BSCCO, 
A||(0) re 1.5 x 1(T 5 cm, s « 1.5 x 1CT 7 cm and A_l/£o « 
10 4 , we estimate 



0.5 



A N (0) 



T 
80K' 



Since [A|| (T)/A||(0)] 2 diverges as T — ► T c we conclude 
that thermal fluctuations may lead to a sizeable reduc- 
tion of the critical current. Note that we study fluctua- 
tions in the linear approximation and, hence, do not take 
into account that the phase perturbations due to ther- 
mal fluctuations may overcome a finite potential barrier. 
The latter process would result in a destruction of a su- 
perconducting current analogous to Josephson junctions 



Observation of smaller values of T c in the stack direc- 
tion, than in direction parallel to the layers, was reported 
in YBCO single crystals with low oxygen content Jl7| . 
Such crystals are expected to have large Ay, which ac- 
cording to (|l5|) is in favor of large fluctuations. However, 
it is difficult to explain by the fluctuation mechanism so 
large differences between critical temperatures observed 
in ref. |jl7j for different directions. 
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